Discussion of a spin-cluster model for the low temperature phase of a'-NaV205 
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We discuss magnetic excitations of a spin-cluster model which has been suggested to describe the 
low temperature phase of a'-NaV205. This model fulfills all symmetry criteria proposed by recent 
x-ray investigations. We find that this model is not able to describe the occurence of two well 
separated magnon lines perpendicular to the ladder direction as observed in INS experiments. We 
suggest further experimental analysis to generally distinguish between models with double refiection 
or inversion symmetry. 



PACS numbers: 75.f0.d, 72.15.N, 7f.20.B, 75.fO.J 

The modeling of the low temperature phase of NaV2 05 
has initiated a discussion of effective spin models for the 
low- lying magnetic excitations. NaV2 05 undergoes a 
phase transition at Tc = 34K ||l[ associated with a lattice 
distortion, charge ordering and the opening of a spin gap. 
In NaV205 the V-ions are arranged in staggered ladders 
along the crystallographic &-axis. While in the high tem- 
peratur phase there is only one V^'^"*" site, recent x-ray 
diffraction studies [||-|| suggest that in the low temper- 
ature phase there are three distinct valence states: on 
every other ladder one finds a zig-zag charge ordering of 
and valence states while on the intermediate 
ladders one finds rungs with two V^'^"*" sites. The struc- 
tural investigation further indicates that the space group 
of the low temperature phase is Fmm2. In the a-h plane 
one finds a doubling of the unit cell along a and b as well 
as mirror planes _La and _Lb. The latter criterion gener- 
ally excludes models with dimcrization along in-line |^] 
or zig-zag chains ||J^ . 

Recent inelastic neutron scattering experiments |^] 
show that there are two close-by magnon excitations with 
a gap of 8.75 meV and 10.65 meV. Both excitations have 
a large dispersion along the &-axis. The magnetic ex- 
change coupling along b has been estimated to range 
between 37.9 meV [D and 60 meV The dispersion 
along a shows only a weak modulation of about 0.5 meV 
which is out of phase for the two well seperated branches. 
Raman scattering experiments observe three excita- 
tions below the two magnon continuum which have been 
interpreted as singlet excitations. Remarkably, the low- 
est excitation has a gap that seems to coincide with the 
gap of the lower branch of the two magnon excitations. 

Based on recent x-ray diffraction experiments in the 
low temperature phase de Boer et al. |^ proposed the for- 
mation of weakly coupled, frustrated spin-clusters along 
the crystallographic &-axis Each cluster (see Fig. 

|I|) contains six Vanadium atoms distributed over three 
ladders and and an overall number of four unpaired elec- 
trons which form a singlet ground state. This spin-cluster 
model is the only proposed model that obeys double re- 
fiection symmetry. 



A previous theoretical study ||f2[ addressed the appli- 
cability of the one-dimensional arrangement to model the 
strong magnon dispersion along the crystallographic h- 
axis. Using a novel cluster-operator theory as well as ex- 
act diagonalization and DMRG calculations the authors 
concluded that there is no parameter regime which would 
reproduce the observed &-axis dispersion. 

In this paper, we present a study of the proposed spin- 
cluster model by means of a strong-coupling expansion. 
We calculate high order results for the magnon disper- 
sions along the 6-axis as well as the leading contribution 
along the a-axis. We discuss several mechanisms to ex- 
plain the occurence of two low 1 ying magnon branches 
without finding a profound supporting argument. We 
point out that symmetries of the magnon dispersions can 
be used to distinguish between classes of models with re- 
fiection and inversion symmetry. Further, we study the 
occurence of singlet states in the spin-cluster model in 
terms of the reported Raman observations. We find that 
there is no evidence for a low lying singlet excitation of 
comparable energy to the lowest triplet excitation. 
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FIG. 1. The spin-cluster model. The filled circles denote 
V*^-ions, the crosses denote V^"'"-ions and the open circles 
denote pairs of V^'^^-ions. 
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The Hamiltonian of the spin-cluster model reads 

H = Ji Si^„ • S2,n + J2 Si „ • S2,n+1 

n n 

+ J' E + ^2,") ■ (S3,n + S4,„) (1) 

n 

n 

+ 84,11 • S'2,n+1 + Si,„ • S'3_„+i) , 

where the Si,„ and S'i_„ denote the four spins on the n'th 
cluster of two neighboring b-axis chains. Ji = (1 + S)J 
and J2 = (1 — (5) J are alternating interactions along the 
6-axis and J3 is the interaction along the a-axis. All inter- 
actions J, J' and J3 are assumed to be antiferromagnetic. 

For an isolated cluster, J2 ^ J3 = 0, we have two 
singlet, three triplet and one quintuplet eigenstates. We 
denote the low-lying eigenstates as follows; 
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(2) 



where and tf^ are singlet and triplet states of the 
spins at sites i and j and V'4 is the = Q component 
of the triplet ■04 . The corresponding energies are Ei = 
-2J' + i Ji, E2 = E3 = -| Ji and E4 = -J' + i Ji. 

The ground state of the spin-cluster is the singlet state 
"01 for antiferromagnetic couplings J and J' with J' > 
i Ji. For smaller values of J' the ground state lies in the 
four- fold degenerate manifold of the states V'2 and -03 • In 
first order, the inter-cluster coupling J2 lowers the energy 
of the tp3 states while leaving the 02 states unchanged. 
One thereby obtaines an effective S = 1 Heisenberg chain 
with a Haldane gap at 7r/(26). 

In the following we will focus on the first parameter 
regime where "01 is the only groundstate. For an isolated 
cluster there are two low energy triplet excitations 03 
and 04. To calculate their dispersion along b we perform 
a strong-coupling expansion around the isolated 

cluster limit treating 



Hi — J2 Si^ra • S2,n+1 



(3) 



as a perturbation. For the moment we neglect the inter- 
chain coupling J3 which we assume to be significantly 
smaller than J due to the longer exchange path along a. 
We have calculated series up to order 7 in J2/J1 where 
the largest system taken into account contains L = 32 
spins. 

The two triplet dispersions are well seperated for large 
values of J' as shown in Fig. |^. Here the obtained se- 
ries are very well converged already in second order. Our 
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FIG. 2. The magnon dispersion for J'=2 J, 5=0.2. The 
solid lines are the result of the 7"^ order strong coupling ex- 
pansion, the open circles are exact diagonalization results and 
the squares denote the results from the cluster operator the- 
ory from |l2[. 



results are consistent with those of the reported exact 
diagonalization and DMRG calculations, whereas the re- 
sults of the linearized Holstein-Primakov approximation 
(LHP) used in the cluster-operator theory |12| disagree 
slightly for momenta in the region around fc 0. 

The perturbative Hamiltonian Hi strongly intermixes 
the two magnon excitations 03 and "04. For small val- 
ues of the cluster exchange J' and the dimerization S we 
find that the two energy bands are very close by for in- 
termediate momenta < fcj, < 7r/(26) as shown in Fig. 



We have calculated effective Hamiltonians for the two 
excitations seperately as well as a combined effective 
Hamiltonian treating the excitations on the same foot- 
ing. We thereby consider the repulsive interaction be- 
tween the two energy bands. While the first approach 
suffers a break-down of the perturbation expansion due to 
small energy denominators in momentum regions where 
the two mixing magnon excitations are nearly degener- 
ate, it allows to assign the nature of the excitation at 
ky = 7r/(26). For example, we find that for J' = 1.1 Ji 
and i5 = 0.5 the lower energy at ky = corresponds to 
a 03-excitation, whereas at ky = 7r/(26) the lower exci- 
tation corresponds to a 04-excitation. In Fig. ^ we have 
used Dlog Fade approximants ||l^ which extrapolate the 
calculated finite series to flatten out the occuring singu- 
larities (dashed lines). 

The second approach allows to explicitely calculate the 
mixing between the two excitation branches. We find 
that there is strong band repulsion for the whole param- 
eter range. For those parameter sets where in the first 
approach the extrapolated dispersions cross we now find 
two well separated branches, but the nature of the exci- 
tation changes depending on the momentum. 

Neutron scattering experiments observe the spin gap at 
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FIG. 3. The ips and tp4 magnon dispersions. In the left col- 
umn the cluster interaction J' is varied for a constant dimer- 
ization 5. In the right column S is varied keeping J' constant. 
A strong mixing of the two branches is observed. Dlog Fade 
approximants (dashed lines) have been used to extrapolate 
the series obtained treating the excitations seperately. The 
solid lines are results of a combined calculation considering 
the band repulsion. For 5 — 0.2 we find only a very limited 
convergence. 



the antiferromagnetic point ky^ = n/b and the zone cen- 
ter ky^ = 0, which would correspond to a i/ja-excitation 
in the spin-cluster model. This assignment enables us to 
determine the ratio of A^ax &t k = 7r/(26) to the gap 
Amin at fc = in order to clarify wether an observed gap 
value of around 10 meV is consistent with an estimate 
of 40 to 60 meV for the spin exchange J. In Fig. 4 we 
show the calculated values at the respective momenta for 
J' — 0.75 Ji as a function of dimerization 6. For small 
values of S the results of the strong coupling expansion 
are less reliable, but it seems that a ratio of 4 does not 
disagree with our results. A tentatively estimated pa- 
rameter set of J' — 0.66 Ji and S = 0.05 would give a 
ratio of around 3 and a spin gap of around 0.25 Ji which 
was doubted in a previous theoretical analysis fl^ . 

Nevertheless, we point out that the shape of the ip^- 
dispersion is rather flat in the vicinity of fcj^ = and ky — 
7r/6 which is due to the strong band repulsion. This seems 
to contradict the experimental observations of a steep 
ascent Further we note that for these parameters 

the minimum of the '04-excitation at ky = 7r/(26) is of 
comparable size to the actual spin gap, but has not been 
reported by neutron scattering experiments. 

The leading contribution to the dispersion of the low- 
lying ^/^a-magnon excitation along a is given by the spin 
exchange mediated by an interaction J3 between neigh- 
boring spin-clusters (see Fig. |l|): 

tab (kx,ky) = eb(ky) 




J' = 0.75 J, 
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FIG. 4. The maximum A,nax of the i/)3-excitation at 
ky — 7r/(2fe) and the gap Amin at fcy = for J' = 0.75 Ji 
versus dimerization 5. The inset shows the ratio Amax/Amin 
for varying values of J' and 5. 



+ 3J.(Jr-L) '°^^^--"^'°^^^--^^- 

The spin-cluster model inherently produces a leading pe- 
riodicity of 27r/a for the dispersion along a. This cor- 
responds to the periodicity observed in INS experiments 
1^,^ and is consistent with a primitive unit cell of area 
2ab as reported by structural x-ray investigations |^ tak- 
ing into account the n/b periodicity along b. The ip^- 
magnon on an isolated cluster consists of a singlet state 
S12 along the 6-axis bond. We therefore expect a strong 
coupling to scattering neutrons at the antiferromagnetic 
point fc^^ = ir/b and only a strongly suppressed signal 
in the vicinity of the zone center at ky"-^ = 0. According 
to (Q) the modulation of the magnon dispersion along a 
is out of phase for these choices of ky. 

Recent neutron scattering experiments |^ report the 
observation of two close by, but well separated, out of 
phase magnon branches. The spin-cluster model at hand 
gives only one low lying magnon branch with the ob- 
served periodicity. This excludes a simple spin-Peierls 
scenario to explain the occurence of two out of phase 
magnon branches, e. g. folding back a single magnon 
branch with double periodicity as it was suggested for a 
zig-zag model in |l^ . 

A way of explaining the occurence of two close 
by magnon branches in this model is to consider an 
anisotropic exchange interaction. An xxz-anisotropy 
could lift the triplet degeneracy in a single branch and 
a doublet branch. Experimentally this splitting of a sin- 
gle triplet branch should result in a 1:2 ratio of neutron 
scattering intensities which has not been reported 
Further, a magnetic field would cause the doublet line 
to split which has not been confirmed by optical spec- 
troscopy JTof . 

The investigation of the spin-cluster model reveals 
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some new aspects that should be covered by future neu- 
tron scattering experiments, namely the experimental ev- 
idence of asymmetric couplings and a verification of the 
basic symmetries found in x-ray scattering. 

Beside a careful investigation of the scattering intensi- 
ties of the two magnon branches, polarized neutron scat- 
tering experiments could give evidence for the splitting 
of a single magnon line, thereby proving the existence of 
anisotropic couplings. 

The verification of the basic symmetries by means of 
neutron scattering experiments is experimentally by far 
less sophisticated. The spin-cluster model is symmetric 
under reflections along the mirror planes _La and _Lb. Ac- 
cordingly, the obtained dispersion eabik^, ky) is symmet- 



ric under transformations of kx 



and kri. 



The experimental data of the dispersion along a for 
ky = and ky — n/b seem to be symmetric under the 
reflection — + — fcx- Nevertheless we point out that 
scanning along an arbitrary value of ky will allow to dis- 
tinguish between models obeying double reflection sym- 
metry and inversion symmetry. For the latter we gener- 
ally expect an unsymmetric dispersion for intermediate 
ky, whereas for models with double reflection symmetry 
we expect a cosine modulation of the amplitude of a sym- 
metric dispersion along a as given in (^. 

For a single spin-cluster, there is one low-lying sin- 
glet excitation 'ip2 which is degenerate with the triplet 
excitation ip^. For vanishing inter-chain coupling J3 the 
Hamiltonian (|l|) is symmetric under a local interchange 
of Ss^n and S4^„. The perturbation operator Hi also con- 
serves this local symmetry. A direct product state of tpi 
states which are even under this symmetry and a V'2 exci- 
tation at some arbitrary cluster which is odd under this 
symmetry couples to a variety of states with the same 
local symmetries, but does not allow the singlet excita- 
tion 'ip2 to move. Though the energy of this excitation 
will be changed, it will not gain any dispersion. As a 
consequence the original degeneracy with the triplet ex- 
citation ip3 is lifted. Neglecting the interaction along the 
a-axis we have calculated the energy of the singlet 1/^2 up 
to order 10 in J2/J1. It turns out that the ratio of the 
obtained energy to the gap of the elementary triplet ex- 
citation is very consistent with the expected value of 1.6 
for the higher singlet excitation at 107 cm"-'^ as measured 
in Raman spectroscopy. 

To explain the occurence of a low lying Raman ex- 
citation at 66 cm-i (8.3 meV) ^ we have estimated 
the binding energies of 2-magnon singlet bound states 
build of the elementary triplet excitations ip^ and ip^. We 
find that in leading order there is no substantial renor- 
malization of these energies and we conclude that these 
bound states are close to the 2-magnon continuum at ^ 
132 cm~^ and are not sufficient to explain the occurence 
of the low lying Raman excitation. Recent ESR stud- 
ies show evidence for a Dzyaloshinskii-Moriya interaction 
below 20 K ||l^,F8| . The inclusion of spin-orbit coupling 



might explain the observation of a Raman excitation de- 
generate to the spin gap. 

In conclusion, we have calculated the magnetic exci- 
tations of a spin-cluster model. We find only partial 
agreement with the experimentally observed spectrum. 
Along the crystallographic 6-axis there are two strongly 
intermixing magnon bands in the parameter regime rele- 
vant for NaV2 05. Due to strong band repulsion we find 
that the lower magnon branch exhibits a weaker disper- 
sion than experimentally observed. Further, the purely 
magnetic model leads to only one magnon branch with 
the observed periodicity along the a-axis. We find one 
low lying singlet state which matches the energy of one 
of the observed Raman excitations. While the prevailing 
models for the low temperature phase of NaV2 05 are in- 
version symmetric, the spin-cluster model obeys double 
reflection symmetry. To generally distinguish between 
these models we propose further experiments scanning 
the magnon dispersions along a for arbitrary values of 
ky. 
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Merit Foundation. 
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